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THE HALF-LIFE OF A DRUG IN RELATION 
TO ITS THERAPEUTIC INDEX 


H. D. LANDAHL 


THE UNIVERSITY OF CHICAGO Toxiciry LABORATORY! 
AND SECTION OF MATHEMATICAL BIOPHYSICS 


A procedure is given for the determination of the effective half-life 
of a drug in vivo with respect to an arbitrary response. If the half-life 
of the drug with respect to toxicity is different from that for therapy, 
the therapeutic index would be expected to depend on the temporal spac- 
ing of the dose. Equations are derived for this variation for a special 
case. 


The stability of a drug in vivo cannot always be measured by 
direct biochemical analysis. It is possible to determine the half-life 
with respect to the production of a particular measureable response 
(mortality, respiratory stimulation, etc.). It is conceivable that the 
half-life with respect to a second response will be different if the sec- 
ond response occurs to some intermediate breakdown product chemi- 
cally different from the one causing the first response. It is also con- 
ceivable that two or more breakdown products will exhibit varying 
potencies toward the first response. In such an instance the determi- 
nation of true half-life is very complex. This latter possibility is 
omitted from the following considerations. What is always measured 
is the disappearance of potency with respect to a particular response. 

Let f(t) be defined as the ratio of potency after time ¢ (in vivo) 
to the initial potency. If doses D, and D, are given ¢ minutes apart 


(D, last) then the effective dose ED is D,f(t) + Dz, where E is meas- 
ured in terms of a reference dose D (e.g., that giving 50% response in 


the population). 
Suppose two identical, and instantaneous doses of magnitude 


(2/3) D are given at ¢ minutes apart. Then (2/3) Df (t) + (2/8)D = 
E(t)D. Since for t = 0, E = 1.33 and fort = 0, H = 0.67 there 
must be some time interval t = T for which E = 1 and hence f(T) = 
1/2. Thus T is the effective half-life. 

1 This work was carried out under contract with the Medical Division of the 
Chemical Warfare Service. 
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The effective dose E produces the given response in a certain per- 
centage of individuals P (or produces a certain percentage of stand- 
ard response in all individuals). A direct plot between P and ¢ can 
be made, and from this plot the half-life T is determined as the value 
for which P = 50 (e.g., 50% mortality) since then E = 1 D (the 
LD-50) and f = 1/2. Percentages of response other than 50%, can 
of course be employed for determination of half-life. 

It is very generally the case that a response to a drug is statisti- 
cally distributed normally with respect to the logarithm of the dose. 
Let o be the standard deviation (here the slope of the log concentra- 
tion-mortality curve) and let D be the LD-50. Then if y is the devi- 
ate of the normal curve, 


o ¥=logi D/D. (1) 


Taking the case in which (2/3)D is Epistles twice at t minutes 
apart, we have: 
(2/3) Df (t) + (2/8)D 


D 


If « and D are known, then f(t) can be determined from 
f(t) = 10 (ov0.176) __ J] . (3) 


If the limiting factor is a monomolecular reaction, and if the doses 
can be considered instantaneous doses then — 


f(t) =2-/", (4) 

and 
o y= log,.(1 + 2/7) — 0.176. (5) 
From expression (5) we can determine for any o the relation between 
y and t/T (the time measured in half-life units) and from this the 
percentage response P as a function of t/T. This transformation has 


been carried out for specific values of « covering the values most com- 
monly observed in toxicological work (Figure 1). 


In this connection it may be pointed out that the slope on linear 
graph paper for the curve relating per cent response to dose in the 
neighborhood of P= 50% is approximately: 


dP/d(logiwD) = 40/c , (6) 

while for f(t) = 2-t/7 
dP/d logit =9/e. (7) 
Because of nonlinearity, the observed slopes are slightly less. Due to 
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FIGURE 1 


Theoretical curves for the per cent response as a function of dose for single 
doses, and as a function of time between application of two doses of magnitude 


2/3 D. The number for each is the value of the standard deviation, «. The 
eurves are on log—probability scale. 


the great difference in the slopes, the standard error of the half-life 
will be much larger than that of the original LD-50. 

In general it is possible to make the experimental determination 
of the half-life and check the LD-50 and o simultaneously. This may 
reveal that the doses D, and D, used were not exactly 2/3 D , in which 
case more general expressions than (2), (8), and (5) are required. 
These are respectively: 


o y= log. [7f(t) + B], (3) 
f(t) =r*10%—s/r, (9) 
77 —logwls +r 27), (10) 


where r and s are the corrected fractions of the LD-50 applied, (7 
being the first). For substances stable in vitro r = s, unless differ- 


ent doses are used. 
- Consider only the case of a monomolecular breakdown so that 


f(t) =2-/7, Let a drug have a LD-50 of D and TD-50 of D’ so that 


its therapeutic index is D/D’. Let the half-life with respect to mor- 
tality be 7 and that with respect to therapy be T’. Suppose that n 
doses are given A units of time apart. Then the effective dose with 


respect to mortality, E,, in units of D, is given by 
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D 
E, —_ se: (1 + 9-4/7 + 2-2h/T a eats pas Ea tiahy 
D 


(11) 


D qd end Q-(n-1)A/T) 
Di—2/7) © 
Similarly the effective dose with respect to therapy, Hr, in units of 
D' is given by 
D 1 = D-(n-1) 4/2" 
Ryee eee (12) 
D' (14 —2“/"") 
Under these conditions the new LD-50 and TD-50, D, and D’, are 
obtained by setting EF; and EL equal to unity. Thus 


— Da—2+7) 
DS; (13) 
(1 Nas Q-(n-1)4/T) 
D'(1—247') 
= —________.. (14) 
(1 =, 2, (eee) 


Denoting by R,, the ratio of the therapeutic index for n doses given 4 
units of time apart to the therapeutic index for single doses (or very 
widely spaced doses), we obtain: 


~—@ — Pa (1 ani D-(n-1) A/T’ ) 
SST er re ee (15) 
4 — 2-4/7’) (1 —Q-m-ayty ” ; 


or for a very large number of repeated doses 
R,= (1—24")/G— 2“). ; . (16) 


For continuous dosage A > 0, and R reaches a maximum or minimum 
value, given by R,, , 


Ra ea > a7) 
Equation (16) is illustrated in Figure 2. 


If instead of equation (4) a more general expression f(t) had 
been used, corresponding expressions could be derived instead of equa- 
tions (11) - (17). The first of those would be of the forms E,D =. 
D[{1 + f(A) + f(2A) +---]. The final result would be modified. If, 
for example, f(t) can be approximated by a linear combination of two- 
exponentials, then the maximum or minimum value of R,, does not 
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necessarily occur at 4 = 0, but at some 4 depending on the values 
of the four time factors. The above results may be applicable to any 
case in which two different responses are initiated by a single agent. 
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on FIGURE 2 
Theoretical curves for the relative therapeutic index (ratio of index for given 
spacing of doses to index for single dose) as a function of the temporal spacing 
of the doses. The time unit is the therapeutic half—life, T’. The curves are on 
log—log scale. 
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b The paper investigates the stability of a metabolizing system sub- 
ject to surface tension forces only. In line with some earlier work of N. 
Rashevsky, the surface tension is considered as a function of the concen- 
tration of the metabolite. Different theoretically interesting cases are 
discussed. The study indicates that surface tension forces alone are not 
likely to produce phenomena similar to those of cell division. 


The work of H. D. Landahl (1942) and N. Rashevsky (1940) 
show the importance of the diffusion forces in cell division. In all the 
cases studied by H. D. Landahl and N. Rashevsky, the cell is assumed 
to have a constant surface tension, which always acts as a stabilizing 
factor, opposing cell division. In his earlier work, however, N. Ra- 
shevsky (1931) considered the possibility that the surface tension 
forces alone may produce a division of metabolizing systems, if the 
surface tension is assumed to be a function of the concentration of the 
metabolites. The lack of proper mathematical technique prevented a 
more detailed study of this mechanism at that time. Subsequent work 
on the role of diffusion forces in cell division threw the question into 
oblivion. However, since the approximation method, introduced later 
(Rashevsky, 1938) enables us now to treat the problem mathematical- 
ly, Professor Rashevsky suggested the present study. While this study 
has now a purely theoretical interest, it contributes to the accepted 
view that surface phenomena alone are not adequate to account for 
cell division. 

Experimental observations on the change of surface tension of a 
solution with changes in the concentration of the solution are well 
known (Bakker, 1928). .For capillary active substances (those sub- 
stances which tend to decrease the surface tension of the solution) 
the change of surface tension with concentration is much greater than 
for capillary inactive substances (those substances which tend to in- 
crease the surface tension). For this reason the former type of sub- 
stance was chosen for this analysis. It is assumed that this substance 
exists in the environment of the metabolizing system at very small 
concentrations and that the relationship between the surface tension 
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of the system and the concentration of the substance can be approxi- 
mated by the linear expression 


y=y,—Ac; fory>0, 
where A > 0. 
With the same notations as those used by N. Rashevsky (1940), 
we have (Rashevsky, 1940, p. 42): 


Since the J, represents the volume forces of the system due to the 
diffusion of the substance in question, and since it is assumed that this 
substance is found in the system only in very small concentrations, 
those volume forces due to its diffusion are negligibly small. This re- 
duces the above equation to: 


1 dr, a 
eet Ls be 3 nV 
The expression for the pressure due to surface tension forces 
(p:) is given (Rashevsky, 1940) by p; = — y(1/R, + 1/R.) where 


FR, and R, are the principle radii of curvature of the surface. Using N. 
Rashevsky’s cell model and notation at the ends of the cell R, = R, = 
72 and at the sides of the cell Ri = 7, and R, = 7.. The pressure at 
the ends due to surface tension is given by p. = — 2y,/7r2 and at the 
sides, ps = —y2(1/r, + 1/72), where y, and y, are the coefficients of 
the surface tension at the ends and sides of the cell, respectively. 

With reference to the Betti equation, (Rashevsky, 1940, chapter 
ili), we have 


1 1 
ey + y¥e= pute (— at = (at sides) 
ry ts 
and 
2 
eLy=— ili si (at ends). 
Therefore, : 


Zz 2 LS | 
I, = — y.7, — S, +4 Y2Ve 3 a5 —) $2, 
To 1," %s 
where S, and S, are the areas of the ends and sides of the cell, respec- 
tively, S, =2a27*, and S:=427,7.. Substituting these values of S, 
and S, in the Betti equation gives 


St 2 CE a 
Sse tara oW A ES Fe ne(—+—)aann |. 
: r 


2 nt 2 
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The volume of the cell, however, is given by (V = 4/327, 7,). The 
above form of J, reduces to 


but a: a = sates I, for surface forces; hence 
Ty dt 3 Vee ‘ ‘ 
1S Ae baal : a 2 
earl ry) ata e 


It must be remembered, however, that y, and y. are given by the re- 
lations y; = 7) — Ac, and y; = yy — AC. where ¢, and c, are the con- 
centrations of the substance at the ends and sides of the cell, respec- 
tively; hence equation (1) becomes 
pees || aco(+++) c 2 | B) 
7, dt 2n a . be 1. (70 ¢) ro 1° 


The next step in the analysis is to obtain c, and c. as functions of 
1,7, D;, De, h, ¢, and q (using N. Rashevsky’s notations), and 
to substitute these functions into equation (2). Using the approxi- 
mation method, N. Rashevsky (1940, chap. i) obtains 

2(DiD. 3 dD;h) © == r,hD Lo 
Cy i 
2D;D-, a 26D;h sr r,hD. 
and 
—-2D,D. + 26Dih + rhDe 
where in the steady state, é=¢, + AqandA=f("1,7,h,6,D;, De). 
For the sake of simplicity, let h > o and let 
A’ =lim aA, ¢ =limé,6“7,, 
how 


oO 


then @ =c, + A’q. Also let c’, = lim ¢, and c’, = lim &. 
h—co how 


Elementary methods yield © 
Bite (2D; + D.) Vie 


fom ese Vert RA ee ae ve) 
ee &9 + 3D. 2(2r.D; > T,D.)Ts Ba (2D, ot B D.) 1°. ( 
or. 7? of. D; -- f. D. 
Cs = Cy + — = - (4) 


SK ef. 
3D. 2(272D; = Dts + (2D; + D.) 12 
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Using equations (3) and (4), expression (2) becomes 


1d ir le 
<= | [mA (9 + x 
Ty dt 2 3 e 


2r.D; + 71D- ) 1 + =) 
2 (275D TDs ie Vee ee z ( Ty Te. 


ate 


e 


(2D; -+D,)%s )] Ze 
2 (2r.D; + 1,De)r, + (2D; + D.)t?n 2) I 8 
Equation (5) is a relation for the relative rate of elongation of 
the cell in the 7, direction due to surface tension forces. When 
the right-hand side of this equation is positive, the cell is elon- 
gating. We shall now determine the conditions which must be imposed 
upon the parameters and variables of the equation, so that elongation 
will occur. 
In order to greatly simplify the mathematical analysis, we as- 
sume that D, = D; = D. This assumption reduces equation (5) to: 


1 dr, ail ee 2r.T*., 2h. 7 
aad Jo A (<2 + > Srreneests SSE wE eam RES. | 
7, db. 2n 3D (4F. > 2n n+ ary 
(= im 1 . 
pee a ) (6) 


Pe fa 3r. 2 
Ml ale Saree] 
3D (4r. + 2r,)7, + 37°, Ns 


At this point, note that when 7, = r, , dr,/dt = 0 , which is what 
we would expect from considerations of symmetry. 


Let it be assumed that the substance is being consumed in the 
cell, ie., g< 0. 


To simplify the algebraic manipulation, consider the following 
notations: 
$= yo— AC 
and 
_ 2r, 1°, \q| A 


C0, SS ese seen nerensiesnenn ner sans FICE 
3D (4r. + 27,) 7, + 87, 
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Substituting w and % into equation (6) gives: 


1 dr, it ee 1 1 2 
a ee ee + 27s * —— SS gl Ra 9 i } =a, ; 
7, at ik Nea epoc o | (- a) Bs aed ole 


1 2 2 


which can be simplified to: 


pial 1 1 ~) + gi 4% 3) (7 
r, at 51 ¢ (= 1» o ( r Ts 2% iE ) 


2 1 


Elongation implies dr,/dt > 0. Therefore since n > 0, we have 


1 1 OYey eu 
Se eae ao 
ete Fk OTs 


Since 7, > 7, , 1/r, — 1/r2 < 0, hence 


which can be simplified to 
(27,—T1) o +o<0. (8) 


CASE I,¢=0. 
Relation (8) becomes 


(27,—T) @ = 0 , 
or, 7, > 272 is a necessary and sufficient condition for elongation. 
CASE II, ¢ > 0. 


Relation (8) becomes 


@ 
(275 —1) gt 1<0, 


or, 7; > 272 is a necessary condition for elongation. 
CASE IlI,¢ <0. 
Relation (8) becomes 
(Ge =r) a 1>0. (9) 


A sufficient condition for elongation is r, 2 2r,. The necessary con- 
ditions are determined as follows: The volume V of the cell is given 
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by V=4/3 a7, r = 4/3 2 1°, where 7) = the radius which the cell 
would have if it were a sphere. Thus 


7?) 
Ty a . 
Po 


ga 


Let M = ——————__ 
3D (yo — ACo) 


Substituting the functions of w and ¢ into relation (9) and letting 
rT, = r,/r*?,, we obtain 


sat (ES ss 2ir*, or Pape ar Ue Mr?, = oon = 0 . (10) 


Relation (10), defined only for 7, > 7 > 7r., is the necessary and suf- 
ficent condition for elongation. In order to investigate the stability 
of the cell when 7, © 7, let r. = (17. — Ar.) where Ar, > 0, then 
neglecting all infinitesimals beyond the first order, relation (10) be- 
comes, after simplifying 


7a Mr, Gs a) 8Ar,) 35 4.575 a) 15Ar, = 0 . 


This implies that for A”, > 0 and Ar, © 0 
4.5 

eee Ee if 
M 


LD os 


(11) 


To summarize, if r, is greater than L, a slight elongation of the 
cell will not result in continued elongation. If, however, ro is below 
the critical value L, a slight elongation will result in continued elon- 
gation. 

Related to the question of “deformation stability” in cells is 
the question of “configuration stability.”” That is to say, for what 
values of the parameters and variables would the surface energy (with 
respect to surface tension) of a cell of given volume be greater than 
the surface energy of the same volume in a two cell configuration? 


This can be determined as follows: 


From considerations of symmetry, when 7, = 72, equation (3) 
is equivalent to expression (4) 


; ~~ ial 
lim c,= lim c’,=c,=c,—_— r2 
T1172 M972 

DisD, Di-D, 


Or 
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Also for a sphere, y, = y2 = yp — Ac, ; therefore, 


»= Alc ala 
read = SS Se 
od) ute 


FE, = Surface energy for a one cell configuration. 
E. = Surface energy for a two cell configuration. 
We require that the two cell configuration be more stable than 


the one cell state. Thus, 


tigate (13) 
but . 


ia Cig 2 |q| 
B,=4aryy=darr| »—A (aso )[, 
and since r for the two cell state is given by r= .87), 
Z . 2 |q| 
Ko—Z An 87 7%, yo —A man rap ° 


Substituting the values of EH, and EF, into relation (138) and sim- 
plifying, we obtain 


gag It 647%) < 0; 
9 D . 0 , 


72, > abate 6<0 
OL ; 
\M| 


7, > =o for ¢>0 
0 GPT . 
[M| 


In summary, when ¢ < 0, the two cell configuration is always 
more stable than the one cell state. When ¢ > 0, however, the two cell 
state is more stable only when 7 is smaller than a certain maximal 
value. By identical methods the case q > 0, ie., production of a sur- 
face active substance has also been studied. 

The following table is a brief summary of all conclusions which 
can be drawn from the analysis: 
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Val f¢ Will cell start Is two cell 
ee to elongate? configuration more stable? 
e=0 Not unless Véa 
1, > 27, 
Not unless 
Substance Not unless perecsri 
Consumed ge ; > 27, r,> # een 
(q< 0) M 
Yes, unless 
0 4.5 Yes 
o< F< 40 
M 
Substance = : 
Produced ¢—0 Ge lark No 
(q > 0) " 2 
Yes, unless Yes 
gees 1, > 2%, 
$< 0 Physically unattainable 


Of the several combinations tabulated above, two are of especial 
interest: choices of g and ¢ which will allow both shape and configura- 
tional stability and choices of g and ¢ which will allow neither. The bio- 
logical significance of these alternatives is as follows: If we imagine 
that surface forces have nothing to do with cell division, then, for the 
best arrangement, the cell division mechanism, whatever it may be, 
should transform the system from one doubly* stable state to another 
doubly stable state. It therefore becomes interesting to know if such 
states are indeed attainable. On the other hand, if we imagine that 
surface forces have everything to do with division we must inquire 
first if doubly unstable cases are at all attainable, it being supposed 
that lengthening is indispensable for cleavage (Gray, 1922). How- 
ever, the possibility of double instability is not the only attribute of 
a satisfactory cell division hypothesis. The additional requirement is 
that the double instability be attained only after the cell has exceeded 
a critical size. Our second problem therefore is to investigate the 
existence of double instability above certain cell sizes. 

What we have called configurational stability is a special case 
of thermodynamic stability wherein the free energy consists only of 


*i.e., with respect to shape and configuration stability. 
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one term, i.e., that for the free energy of the surface. It can be seen 
from Table 1 that for g < 0, the system is thermodynamically stable 


only when ¢ = 0 and 0 < ~» = V—2.34/M . It may be asked whether 
there are various shapes (i.e. various equilibrium points with respect 
to changing shape) consistent with this restriction. Equilibrium with 
respect to changing shape requires that relation (10) be an equation. 
Putting § = r./r, into relation (10) as an equation, we find 
( Sit Bi4~2 
Tr) = | ———-_ 
2M f? (1 — 2 6?) 

The general shape of this curve can be deduced from the fact that, 


becomes infinitely great at the left end of the curve (f > W¥ and 
To =x > 0 when 6 > 1. At these points, the derivative 


dry 1 
ap -2r, BIL — 2B) 


1/2 ' 
) mvs — Bol. (14) 


38 (86° +2) —2r,M (1—5p%) (15) 


is respectively large and negative ($ > W4) and negative when f > 1. 
By setting expression (15) equal to zero, we obtain 
f(f) =3 2? —14  — 12. B'+2=—0. (16) 

The quantity f(f) is an odd function and has two changes in sign, 
hence it can have but one zero in the range of (6. This means that 7 
plotted against £ is roughly of the shape of a half U having an asymp- 
tote at 6 = ¥/4 and a positive value « at B=1. 

The coexistence of shape and thermodynamic stabilities may now 
be ascertained by finding if it is possible for the curve (14) to inter- 
sect the straight line parallel to the f—axis, 


Zoku 
To = ( = 71 ) : (17) 
By equating expressions (14) and (17) we obtain 
3 6° — 9.36 f° + 4 B* + 4.68 6 + 2=—0 (18) 


as the condition for the intersection of expressions (14) and (17). 
It can readily be seen that equation (18) has no positive roots for 
WH < B <1, hence the shape, and the thermodynamic stability can 
never coexist for this case. This result need not preclude the occur- 
rence of thermodynamically unstable shapes (non-spherical) if the 
process of division, whatever be its nature, has a high energy of ac- 


tivation. 
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For the case when g > 0 and ¢ = 0, the cell is in equilibrium with 
respect to shape change when relation (8) becomes an equation with 
¢ = 0, 1e., 

p= (4). (19) 


Since for this case, the system is always thermodynamically 
stable, all cells satisfying equation (19) will be stable with respect to 
both shape and configuration. 

The first question posed above can therefore be answered in the 
affirmative. Double stability throughout a range of shapes is un- 
equivocally attainable for gq > 0 and ¢ = 0. There is a slight chance 
that it would be attainable for gq < 0;0 < ™ < VY—2.34/M . 

The answer to the second question is obvious from inspection of 
Table 1. The only possible combination for double instability (q > 0; 
¢ > 0) operates below a critical state, rather than above, as would be 
required of a satisfactory division mechanism. 

Conclusions. It is generally agreed upon that any analysis of 
the mechanism of division must take into account several contribut- 
ing factors, whose resultant effect is the phenomenon of cell division. 
Observation on the structure and division of different “kinds” of cells 
leave no doubt that those various contributing factors are rarely if 
ever present in the same “proportions” in any two kinds of dividing 
cells. 

A question early raised in connection with this work is: which 
of these factors, if any, is the dominant one in cell division? One ap- 
proach to the answer of this question is the analysis of each of the ap- 
parently important factors separately with the assumption that it 
alone could at least qualitatively account for the phenomenon. In the 
present paper, using the above approach, we start out with the as- 
sumption that surface tension is this dominant force. It has been re- 
quired simply that surface tension effects would account for the elon- 
gation and thermodynamic instability of cells above certain critical 
sizes. Any additional mechanisms which would be required to explain 
the transition from the thermodynamically unstable elongated cell 
to the two cell state have been tacitly assumed to operate. The analysis 
leads to the conclusion that surface tension forces, in spite of these 
additional sympathetic assumptions, do not account for the process 
of cell division. This conclusion is in complete agreement with other 
theoretical and experimental evidence, chief of which is the extremely 
low tension at cell surfaces measured by recent workers (Bourne, 
1942). 


It should be noted, in closing, that the conclusions to which we 
have been led are necessarily dependent upon the nature of the ideal 
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system we have assumed. A subsequent analysis should certainly pro- 
ceed without the following of the assumptions made here: infinite per- 
meability to the surface active solute, zero adsorption of the solute 
at the cell boundary, and constant rate of metabolism. Ultimately, of 
course, the theoretical scheme must be patterned after what is known 
regarding the metabolism of surface active agents within cells. 

The many valuable suggestions offered by Dr. M. Morales and 
Professor N. Rashevsky greatly facilitated the writing of this paper. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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Evidence is presented to show that despite their effectiveness in 
producing motion sickness the actual frequencies and accelerations in- 
volved in ship motion are rather low. It is suggested that the disparity 
is due to the fact that ship motion is such as to cause an asynchrony 


between signals from the labyrinth and those from other motion recep- 
tors. 


Introduction. In order to obtain an incidence of motion sickness 
among the test population comparable to that aboard ships, experi- 
mental devices operate at peak accelerations of 1g (980 cms sec”) or 
more, and at frequencies in considerable excess of ship oscillations. 
A comparison of these values with real ship motions can be made 
assuming that the latter are simple harmonic in both roll and pitch 
(Fig. 1). Such computations have been made for certain positions 


ROLL PITCH 
(2) 
SECONDS 
Ot ae Soe eed ee ee a SS eet 
O 4 8 [2 16 20 24 28 32 O 4 8 l2 16 20 
FIGURE 1 


Simultaneous records of roll and pitch taken from unpublished observations 
aboard the U.S.S. Washington (BB 56). 


aboard several classes of naval vessels (Table 1). The amplitudes 
used are from sallying data, and can be considered as absolute upper 
limits; nevertheless it will be observed that the average peak accel- 
erations developed are well below 1g. This disparity in the effective- 

i i ded to give the impression that the peak value of the accelera- 
tion ester en rieal Paks only that it is not the sole factor. The order of 


in Table 1 (descending order of accelerations) is almost precisely the or- 
ieee a ee ee ee in producing seasickness among the crews. 
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TABLE 1 
| Period Amplitude Period Amplitude ‘ 
Wert Te | (SERS | Seren | oem) | Begs [deere 
Destroyer 9.5 25 5 5 .26 ¢ 
LCI 4.5 25 3 5 WALES 
LST 5.0 15 4 3 19¢ 
Destroyer Escort 8.0 25 5 5 14¢ 
Heavy Cruiser 12.0 20 7 4 14¢ 
Carrier 16.0 15 7 4 13.¢ 
Light Cruiser 12.0 20 Hj 4 silane 
Escort Carrier 13.0 20 ff 4 12¢ 
Battleship 14.0 12 use 6 3 12¢ 
Transport 14.0 25 7 4 .08 ¢ 


This table is constructed entirely from unrestricted information. The periods 
and amplitudes are average values kindly furnished the author by the Bureau of 
Ships, Navy Department. The “average acceleration” is the arithmetic mean of 
the roll and pitch contributions at three points aboard the vessel: at the bow on 
the water line, abeam of midships on the water line, and high on the foremast. 


ness of experimental devices and real ships naturally suggests that 
the shipboard motion has certain additional features not yet incor- 
porated into the devices, or, at least not in the proper quantitative 
manner. Ultimately such features must be related to the characteris- 
tics of the peripheral and central mechanisms by which the organism 
perceives the motion. This paper is concerned with features which 
appear to be directly related to the mechanics of the labyrinthine re- 
ceptors. It must of course be recognized that in addition to the ques- 
tion of effective stimulation, motion sickness poses the problem of 
how the stimulation operates to fire off the visceral responses. This 
second question undoubtedly involves central phenomena (E. A. Spie- 
gel, 1944), and no attempt is made here to establish a liaison between 
the two questions. 

There can be little doubt, especially since the study of A. Sjéberg, 
(1931), that labyrinthine, visual, and proprioceptive stimulation each 
play partial roles in the etiology of seasickness. A. Sjéberg demon- 
strated this by excluding these sources of information, one at a time, 
but it is reasonable to suppose that asynchrony of one receptor with 
respect to another would likewise affect adversely the ability to resist 
the sickness. The asynchrony might be described as a situation in 
which “coordination is difficult” and although this description is vague 
with respect to the central mechanism of “coordination,” it expresses 
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incontestable experimental facts. The asynchrony proposed here as 
an etiological factor is chiefly that which can arise between labyrin- 
thine mechanisms and those of other receptors. More specifically, there 
1s comparatively little time lag between the application of a variable 
external stimulus and the variable physiological response in the sen- 
sory fibers from the eye or from proprioceptors. If 6 is the instanta- 
neous angle of a periodic roll or pitch, it may be assumed that the 
frequencies in such fibers are in phase with 6. This may or may not 
be the case for the labyrinth, as will be discussed below. 

Mechanical Considerations Regarding the Canals. The discus- 
sion of this section for the most part derives its factual information 
from the beautiful experiments of G. Schmaltz (1931) on man, W. 
Steinhausen (1933) on the pike, and O. Léwenstein and A. Sand 
(1940) on the ray; therefore at the outset the dangers attending trans- 
specific analogy must be acknowledged. 

G. Schmaltz and W. Steinhausen have established that for me- 

chanical purposes a semicircular canal—in particular, the horizontal— 
can be treated as a rigid toroidal container filled with a fluid of def- 
inite viscosity. It is here assumed that the remaining canals can like- 
wise be schematized.? Into the lumen of the canal, at the ampulla, ex- 
tends the cupula, much as a gate hinged at the crista, containing the 
nerve endings and with them constituting the end organ. 

Consider now the rotation of the canal about an arbitrary axis 
normal to its plane. Two simple but important facts can be readily 
shown: (a) The forces acting to rotate the fluid within the canal are 
determined exclusively by the components of the acceleration normal 
to the radius of gyration, and (b) The effect of these forces on the 
fluid will be the same as though the canal were being rotated with a 
similar motion about its own axis of symmetry. These circumstances 
make it easy to describe the motion of the endolymph as the whole 
organism is being rotated about a distant axis (the usual experimental 
situation). 

If w and w, are respectively the angular velocities of the canal 
and of the contained endolymph, and uw is the coefficient of friction 
calculated from Poiseuille’s Law, then it can be shown that, 


t= u(@— or). (1) 


If the viscosity coefficient and density of endolymph are taken as equal 


2 For the thesis developed in this paper it is only necessary that the pattern 
of discharge in all three canals be affected in the same way by the physical 
stimulus. However, it should be emphasized that corresponding central effects 
may be (and, on physiological evidence, likely are) different for the three canals. 
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to those of water, and the radius of the canal is taken as 2 X 10° cms, 
then uw is 200 sect. 

Equation (1) can be integrated for various motions imposed on 
the canal to discover the corresponding motion of the endolymph. 
Three of these are of experimental interest: I. The canal is rotated ata 
constant angular velocity, starting from rest. Il. The canal is acceler- 
ated angularly at a constant rate, starting from rest, and III. The 
canal is oscillated in angular simple harmonic motion (6 = 4 sin wt). 
In each of these cases there will be a certain pressure over any cross 
section of the endolymph, i.e., the pressure which would act on any 
structure fixed rigidly to the canal and extending into the endolymph. 
These pressures, P;, Py, and Py, are given by the product of P,, a 
constant common to all cases, a constant determined by the particular 
experimental conditions, and a time dependent factor, thus, 


P= Pf, fer 3 (2) 
Py=P,P.(1—e#'), (3) 
ll 
Py=— P,P; cos(wt+¢), (@¢=tan'?—). (4) 
@ 


Recalling that wu is 200, it is evident that in case I, the pressure on any 
structure will vanish very fast. In case II it will attain a constant 
value very fast. The angular velocity, w , imposed by ship motion is 
of the order of 1 sec“ so that in case III ¢ is an angle very near 90°, or 
Py, = P,P; sin at, ie., the pressure is almost exactly in phase with 
the imposed motion. The third case has never been tested experi- 
mentally, but the conclusions regarding the first two have been verified 
by G. Schmaltz. W. Steinhausen recognized that this general syn- 
chrony between the imposed motion and the endolymphatic motion 
precluded the explanation of persistence phenomena, such as a 20 
second nystagmus, in terms of fluid movement. In resolving the dif- 
ficulty he suggested that the site of persistence phenomena might be 
the cupula itself, in the sense that it deflected as a damped vibrator. 
W. Steinhausen convincingly proved this hypothesis by direct obser- 
vation on the pike, showing that under experimental condition I the 
cupula is rapidly deflected from rest and returns slowly to the origi- 
nal position some 20 seconds later. Unfortunately Steinhausen’s 
analysis appears to assume a “rigid linkage” between endolymph and 
cupula, and thus to contradict his experiments. Rectification is easy, 
however. We need only consider the deflection of a damped vibrator 
under the influence of the imposed force functions given by equations 
(2) — (4). Strictly speaking, the force on the cupula is proportional 
to these expressions and to the cosine of the angle through which the 
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structure is deflected, but if we limit ourselves to deflections of less 
than 25°, the cosine can be taken as unity. Integration of the equa- 
tion of motion then leads to expressions for cupular angle of deflection 
p(t) for each of the three experimental situations. Since there is no 
available data on y(t) itself, but only O. Lowenstein and A. Sand’s 
measurements on the change in resting frequency of discharge in the 
ampullar nerve, »(€) — », the theoretical predictions cannot be test- 
ed without an additional hypothesis regarding the relationship between 
cupular deflection and nerve response. The assumption made here and 
justified a posteriori is that the per cent change in frequency of nerve 
discharge is proportional to the angle of deflection from rest. If this 
is the case, then »(t) and y(t) should be curves of the same general 
form and they appear to be so (see Lowenstein and Sand’s discus- 
sion). The analytic relationship would be »(t) — ») = h» w(t), where 
h is a proportionality constant. 

Taking 24 and x’ as gross coefficients of friction and elasticity re- 
spectively, it can be shown by straightforward methods that, 


P, jee u 
a (en? — ent) , (5) 
Ke 2A bP Ty 2 
where 7; =—A + V/22— 2 and 7, =—A — Vi? — k?. This function has 


the well known properties of a rapid ascent, passage through a maxi- 
mum, and an elongated return to zero, as demanded by the experi- 
mental data for the case. When fitted in the usual way to O. Lowen- 


stein and A. Sand’s data (p. 270, Fig. 9), it is found that 7, = — .058, 
and r, = — 10.6, and these values are independent of the maximum yw 
attained. Using this 7, and 7, it is found next that, 

iy = 1.78 P, Po (i — e-°') = 0944 P, Pt, (6) 


the linearity holding very accurately for the first five seconds. Since 
P, = a/u, equation (4) states that within the first five seconds y is 
proportional to the angular acceleration and to the time. Again this 
is in good agreement (ibid p. 269, Fig. 8). However, curve fitting is 
not weighty evidence per se. A good cross check can be obtained by 
taking P, from the data of case II and using it to calculate the maxi- 
mum y in case I. This calculated maximum (28.4 sec’) is in good 
agreement with that given by the experimental data of case I (26.0 
sect). 

With these damping characteristics of the cupula extracted from 
the data for cases I and I, we can now predict the pattern of the dis- 
charge in the ampullar nerve when the canal is subjected to the 
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sinusoidal motion, 6= 6, sin wt. For case III the equation of motion 
gives, 


(vy — ») =Constant X cos (wt + ¢ + &), (7) 
where 
Lak 
= tant ———. (8) 
Ge Kb 


It has already been stated that ¢ = 90°, consequently if € were zero, 
equation (7) would predict that the ampullar discharge would be in 
phase with the angular displacement, @. On the other hand, the dis- 
charge would be maximally out of phase for § = 90°. Thus € meas- 
ures the asynchrony. Equation (8) shows how & depends on the elas- 
ticity and friction of the cupula and (noting that the period of the im- 
posed vibration, 7 =2z/w) on the period of the imposed vibration. Fig- 
ure 2, which is a plot of & vs. T using the experimental values of r, and 
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FIGURE 2 


Plot of the period of the imposed motion vs. the phase angle, £ of the cupula 
§ = 90° indicates maximum asynchrony. The ranges of roll and pitch for major 
naval vessels are indicated by arrows. 


r, and the relations 7, + 7, =— 22; ryr. = x’, shows that asynchrony is 
maximal for a range of periods precisely equal to the periods of roll 
and pitch of large vessels. When T is increased the asynchrony drops 
too, although slowly. (However when T is increased the amplitude is 
decreasing so that very likely there is no problem anyway.) Thus it 
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appears, if the present suggestion is correct, that so far as the canals 
are concerned, vessels are constructed so as to maximize the confu- 
sion of their passengers. 

Mechanical Considerations Regarding Component Receptors. The 
discussion of the first section has been concerned exclusively with the 
mechanics of semicircular canals. The remaining labyrinthine recep- 
tors are presumably the utricles. Their precise mechanism appears 
to be unknown, but their anatomy suggests that, unlike the canals, 
they are stimulated by one of the principal components (horizontal 
and vertical) of the imposed motion. It will be shown in this section 
that these components are in phase with @ for certain amplitudes and 
positions aboard the vessel and are out of phase for others. 

Consider the roll of the vessel schematized in Figure 3, where P 


FIGURE 3 
Motion of a Point, P, which is at position angle, y, when the hull is level. 


is some point fixed on a cross section through the hull. F is the dis- 
tance from P to the axis of rotation, O. When the vessel is at the 
equilibrium position, the position angle of P is y. As the vessel rolls, 
therefore, the coordinates of P are given by: 


x—=Reos(é+y), (9) 
= Rsin(¢ + y). (10) 


The horizontal and vertical components of the force are of course pro- 
portional to d?x/dt? and d*y/dt*, and these can be found from equa- 
tions (9) and (10) under the assumption that the ship rolls in angu- 
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lar simple harmonic motion, 6 = 6 sin wt. The times at which these 
force components reach a maximum or minimum can be found by set- 
ting equal to zero the third derivatives of x and y, and solving for 
the corresponding t. When this is done, it is found that for any one 
position angle, the horizontal force stays in phase with @ unless the 
amplitude, 6), exceeds a certain critical value 6*. The critical ampli- 
tude varies with y, i.e., with position aboard the vessel. Likewise, the 
vertical force stays in phase with 6 unless 6) exceeds 6**, in which case 
it falls out of phase during the first half of the cycle. If 4 is increased 
still more, to the point where it exceeds a second critical value 6°’, 
then an asynchrony develops during the second half of the cycle as 
well. Again, 6** and 6*** vary with position aboard the vessel, i.e., 
with y. These relationships are best understood from Figures 4 and 
5. Illustrative wave forms are shown in Figure 6. 


40 


IN THIS ZONE THERE IS AN INTERMEDIATE 
MAXIMUM ON THE NEGATIVE SIDE OF THE . 
WAVE. 


30 


20 


IN THIS ZONE THERE IS NO INTER- 
MEDIATE MAXIMUM ON EITHER SIDE 
OF THE WAVE. 


CRITICAL AMPLITUDE, IN DEGREES 


ie) 10 20 30 40 50 60 70 80 90 
POSITION ANGLE, % IN DEGREES 


FIGURE 4 


Critical angle (9*) for fundamental changes in the horizontal acceleration 
wave form, as a function of position angle, (y). 


It appears from the foregoing considerations that a receptor sen- 
sitive to horizontal or vertical components of the force developed in 
a roll may or may not be synchronous with receptors such as the eye 
or the proprioceptors, which are almost exactly in phase with @. In 
this respect such receptors resemble the semicircular canals. An addi- 
tional complication arises, however, in that the asynchrony for com- 
ponent receptors varies with position aboard the ship. 

Discussion. If the hypothesis made in this paper to the effect that 
asynchrony between various receptors is one of the major factors 
which induce seasickness, then it appears that shipboard motion may 
owe some of its effectiveness to the fact that it imposes an asynchrony 
both on receptors which respond to circular flow of fluid and to those 
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Critical angles for fundamental changes in the vertical acceleration wave 
forms, as a function of position angle, y. 
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FIGURE 6 
Schematic theoretical time variations of linear acceleration to illustrate wave 
forms for different combinations of amplitude and position angle. Upper curves 
are for horizontal components; lower curves for vertical components. The verti- 
cal dimensions have all been set so as to keep maximal heights the same regard- 
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which respond just to particular components of the imposed motion. 
In the former case it appears further that, with pitch—roll periods as 
they are and mechanical constants of the cupula as they are, ship- 
board motion actually maximizes the asynchrony. 

In the past, rocking devices, swings, and elevators designed to 
study seasickness experimentally have seldom been constructed so as 
to reproduce accurately real shipboard motion. Accurate reproduc- 
tion may be rather essential in view of the foregoing results. Records 
of ship motion (e.g., Fig. 1) can be made routinely aboard naval ves- 
sels, using fire-control equipment? and these wave forms may then 
be reproduced by the device. Aside from obtaining wave forms which 
cause asynchrony, there are preserved other features which may be 
of etiological importance. Among these is the fact noted in unpub- 
lished measurements of the author that a stormy sea as contrasted 
with a calm sea can be characterized just as well by a greater spread 
in amplitudes as by a greater average amplitude. This erraticity 
would seem to hinder any progressive adaptation. 

In closing may be mentioned a problem suggested by cupular me- 
chanics but one which may have some importance in general sensory 
physiology. The integral, S|y|dt, measures in a sense the total stimu- 
lation which the canal receives, and it may be asked, what wave form 
of the wmposed force will maximize this integral? This is clearly a 
problem in the calculus of variations, but it is one of whose solution 
the writer is not aware (and one which the writer has not solved). 
The solution must be something akin to sinusoidal motion for the rea- 
sons which O. Lowenstein and A. Sand discuss. A related problem is 
to find the wave form of the imposed forces such that py will have a 
given form. This can in many cases be done, and the results may be 
of some use in further experimental studies. 

In the preparation of this manuscript the writer has had the 
helpful counsel of Professor Spiegel of Temple University and Pro- 
fessors Andrews and Kleitman of the University of Chicago. Certain 
of the work in the second section appeared in a Navy report with Mr. 
James Birren (1945). 
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An integral equation is established applicable to many types of me- 
tabolizing systems. The equation is applied to first order chemical re- 
actions and to some biological systems. 


A metabolizing system usually presents so many variables that 
a mathematical description is either impossible or crudely approxi- 
mate. By neglecting some of the details of such systems, a satisfac- 
tory description, preserving important features, is possible through 
integral equations. The types of equations encountered are fortunate- 
ly well known in mathematics. Thus if the description is adequate 
for the biological problem, the treatment will open a broad range of 
mathematics to the biophysicist for treating metabolizing systems. 

The Integral Equation: For our purposes, a metabolizing system 
is any localized region where a substance or material is being con- 
sumed, produced, transported, modified, or stored. The system 
may be performing all or several of these functions simultane- 
ously. In general, the system must have some defined boundary in 
order that experimental results may be obtained. In theory, however, 
a physical boundary need not exist. . 

A general metabolizing system may be pictured as an irregularly 
shaped region of complex structure into which material is entering, 
in which some material is being produced or consumed or stored, and 
some may be leaving by diffusion through the walls or by flow through 
openings. We introduce the following symbols: 


M (0) =amount of the metabolite present initially, t = 0. 


M(t) = amount of the metabolite present at time, t. 
R(6@) =rate at which the metabolite is accumulating in the sys- 


tem. 
F(t) =the “metabolizing” function, to be defined presently. 


In this region of irregular geometric shape, we are interested in 
some metabolite originally present in amount, M(0). This metabolite 
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is accumulating at a rate R(@) so that the amount accumulated in 
time interval dé is R(0)d6. This additional amount suffers the fate 
of the original; it is metabolized. We introduce a “metabolizing” 
function defined as that function which will multiply the original 
amount, M(0), to yield what remains of it at time, t. That is, the 
amount remaining from the original amount at time, t, is M(0) F(t) ; 
the amount entering in the interval, 6 to @ + dé, which remains at 
time, t, is F(t) R(6) dé. The total amount present at time, ¢, 
from both these contributions is M(t) which is seen to equal 


M(t) =M(0) F(t) + [ R@ Fe—#) dé. (1) 


This equation is similar to some treated in population studies (V. 
Volterra and U. D’Ancona, 1935; A. J. Lotka, 1939) and is almost 
identical with one used by R. V. Churchill (1944) in another problem. 

Whether equation (1) is an integral equation of the first or sec- 
ond kind depends upon the information we have of the system. The 
most likely situation experimentally is that we know M(0) and can 
determine M(t), then if we can control or determine R , the equation 
is a Volterra integral equation of the second kind in F(t). Occa- 
sionally the information will be such that R is the unknown, then 
equation (1) is a Volterra integral equation of the first kind in R (H. 
Margenau and G. M. Murphy, 1943). Finally if M(0), R, and F are 
known, the equation can be integrated either directly or numerically 
to find M(t). 

Applications to some chemical and biological systems: A large 
number of biological and chemical systems are characterized by 
the fact that the rate of reaction of the substance of interest is 
directly proportional to the quantity of the reacting substance 
present. This type of reaction is designated in chemical kinetics 
as a first order reaction. As a numerical example, we can consider the 
decomposition of dibromosuccinic acid in hot water investigated by 
Van’t Hoff and treated by F. L. Hitchcock and C. S. Robinson (1936). 

All first order reactions are very easily treated by equation (1) 
with the following designations: 


M (0) =a known constant. 
F(t) =F) =1. 
R(0) =—k M (6). 


Equation (1) becomes 


M(t) =M(0) —k fw do. (2) 
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Equation (2) is of the type 
b 
(a) =f(e) +4 | K(@,2) 6(2) ae. 


This Volterra equation of the second kind is solved by bt Liouville— 
Neumann series (Margenau and Murphy, 1943), 


o(«) = > Aha (2%) 


n=0 


provided 
do(%) =f (x); 4,(x) = | K (% 2) do (Z) dz; -++gn (4) = 


HM BAD 2) tas ee 


In our example, K(x ,z) =1,i=-—k, and ¢.(x) = M(0). 
Introducing these expressions and integrating, we have 


te 


M(t) =M(0) > 


n=0 


=M(0)e*. (3) 


Equation (3) is the usual exponential expression which is more easily 
obtained in this instance by stating the relation in differential equa- 
tion form and integrating. With k = 0.030, the equation is an ade- 
quate description of the decomposition of dibromosuccinic acid in the 
Van’t Hoff experiment. 

First order reactions in certain systems, however, seem possible 
of treatment in much more elegant manner through equation (1) than 
through differential equations. As an example, let us consider a sys- 
tem in which the reaction is of the first order and at the same time 
the material is flowing out constantly such that the “metabolizing” 
function is of the form F(t) = (1 — ct). Using the designations as 
in equation (2) we have 


M (t) =M (0) (1—ct) —k iO [1—e(t—6)] dé. (4) 


This integral equation may be solved by the Liouville-Neumann series 
as used in equation (2), however, it is desirable to introduce the pow- 
erful mathematical technique of the Laplace transfornt which is emi- 
nently suited for treating many of the problems described through 
equation (1). 

The integral on the right of equation (1) is a special type re- 
ferred to as a “Faltung” integral (literally: a folding integral) by 
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G. Doetsch (1943). Equation (1) may be written. 
M(t) =M(0) F(t) + FR, 


where F°R is the abbreviation for the integral. If we multiply through 
by exp (—st) and integrate from 0 to «, we have 


fue) e-*t dt =M(0) xo e-8t dt + 
eis (5) 
{ f R(0)F(t—6) et dt do. 


The conditions on the functions in order that these integrations are 
permissible are discussed by G. Doetsch (1943). We shall assume 
that the integrations are allowed. By definition the integral on the 
left of equation (5) represents the Laplace transform of M (¢t) 
which we shall call m(s) ; the first integral on the right is the trans- 
form of F(t) multiplied by M(0), that is, M(0) f(s). The double 
integral on the right is the transform of the “Faltung” integral. G. 
Doetsch (1943) shows that the transform of F*R is equal to the prod- 
uct of the individual transforms, or 


L{F*R} = L(F) -L(R) =f (s) r(s). 
Hence equation (5) may be written as 
m(s) =M(0) f(s) + f(s) r(s), (6) 


where the lower case letters represent the transforms of the corre- 
sponding functions in capitals. An alternate statement is that equa- 
tion (1) transforms into equation (6). If we seek f(s) we have 
simply 

m(s) 


M(0) + r(s)° 


It may happen that the particular form of the right side of equation 
(7) is one of the many Laplace transforms which have been tabulated 
or that it may be expressed in standard form through algebraic manip- 
ulation. If either circumstance exists, we may read directly from 
tables of transforms given, for example, in G. Doetsch (1948), H. S. 
Carslaw and J. C. Jaeger (1941), R. V. Churchill (1944), or N. W. 
MacLachlan (1939), the function F (t) which we seek. If the right 
side is not a standard form, integration in the complex plane is re- 
quired. The inverse transform is 


f(s) = (7) 


1 
Os ser {ec Hs) ds. 
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Details of the integration may be found in the texts mentioned above. 
. Applying the Laplace transform technique to equation (4), we 
ave 


fier M (t) dt 


= M (0) { e*t(1— ct) dt —k ffm (6) [l1—c (t—6)] e-** do dt. 


These integrations are easily performed or the transforms may be 
read directly from a table: 


m(0) em(0) bape 
ee eae 


m(s) = : 
Ss Ss 


This equation may be written in terms of m(s) with the right side 
in a standard form, 


A B 
m(s) =m(0) eee 
S10 a 
where 
a=i(k+ Ve +40), ieee 
Za—k 
see be dia 
b=4(—kti VP+9kc), B= 
2a—k 


Consulting G. Doetsch’s (1943) table again we find the inverse trans- 


forms to be 
M(t) =M(0) [Ae + Be], (8) 


Equation (8) describes behavior of a system in which a first order 
reaction is taking place and, at the same time, the reacting material 
is being transported away at a constant rate. It is evident, too, that 
equation (8) reduces to equation (3) whenc=—0O. 

The system described by equation (8) seems useful in biology. 
For example, the breakdown of glucose in tissue may be considered 
with some approximation to follow a first order law. Hence a piece 
of muscle tissue in vitro bathed continuously in a glucose—free solu- 
tion would be an example of such a system where M(t) would refer 
to the concentration of glucose. 

A large number of biological systems under normal conditions are 
characterized by the fact that the amounts of certain metabolites re- 
main constant. In terms of equation (1), we have for each such me- 


tabolite 


164 METABOLIZING SYSTEMS 


MUL— F(t] = {| RO) F—#) 46, (9) 


where M is the constant amount present. 

The procedure from this point is determined by what additional 
information we have. For example, we may know or suspect that the 
rate at which material enters the system is constant; our problem be- 
comes one of determining the F-function. The resulting equation is 
easily shown to have the solution F(t) =exp (—R/M)t. Thus a sys- 
tem in which the amount is constant is characterized by an exponen- 
tially decreasing metabolizing function. The converse of this prob- 
lem is also of biological interest. It is interesting to note that a me- 
tabolizing system with a constant amount and a constant F—function 
is impossible, for assuming that R(t) is a well-behaved function, this 
system would require 


t 
i R(6) dé=constant, 
0 


which is not possible. 


‘The rate of change of the amount of the metabolite may be found 
from equation (1) by differentiating, 


aes oF (t—é 
pee mos fixe (¢—8) 


a 


ont (hE) Pe 


bse tele, is of additional interest, ae it may be used as the inte- 
gral equation of the second kind for determining R(t). 
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=~ ERRATA 


x 2 In the paper “Chain Processes and Their Biophysical Appli- 
cations: Part I. General Theory” by I. Opatowski, Volume 7, pp. 
et on page 177 equation (38) should read: 


4 Vint k= —- Yoj-2" ae , 
(40) should read: 3 
| Yoel (ita eres: oe = Ye Yin Yrs” ce oe Y ne 
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1. The Bulletin is devoted to publications of research in Mathe- 
matical Biophysics, as contributing to the physicomathematical foun- 
dations of biology in their most general scope. 

2. Papers published in the Bulletin cover physicomathematical 
theories as well as any other mathematical treatments of biological 
phenomena, with the exception of purely statistical studies. 

3. Mathematical studies in physics or in borderline fields in 
which a direct connection with biological problems is pointed out are 
also accepted. 

4, Emphasis is put upon the mathematical developments, but a 
description and discussion of experimental work falls also within the 
scope of the Bulletin provided that description or discussion is made 
in close connection with mathematical developments contained in the 
same paper. 

5. Outside of the scope of the journal are papers of purely sta- 
tistical nature or papers concerned only with empirical equations. 
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right consecutively. Do not use such notations as “equation 2a” or, 
“equation 5” etc. References should be all given at the end of the 
paper, arranged in alphabetic order by authors and for each author 
chronologically, following exactly the style used in the Bulletin. In 
text, reference should be made by giving in parentheses the name of . 
the author followed by the year of publication. In case of several pub- 
lications by the same author in the same year, use notations “1940a”, 
“1940b’’, etc. 

In writing equations, slanted lines should be used wherever pos- 
sible. 

Every paper is to be preceded by a short abstract. 

Drawings should be prepared in a professional manner on white 
paper or tracing cloth following as closely as possible the style used 
in the Bulletin. They should be drawn approximately twice the scale 
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the Editorial Office according to author’s instructions, at cost. 
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